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WSTĘP 


Plechy glonów z rodziny ramienic są zbudowane z dużych komórczaków 
— międzywęźli, które oddzielone są od siebie pojedynczymi warstwami 
mniejszych komórek — węzłami. Przypisuje się im zwykle funkcję łącz- 
nika, zapewniającego transport substancji, zwłaszcza substancji mineral- 
nych i produktów fotosyntezy od jednego międzywężźla do sąsiedniego [6, 7]. 

Gdyby jednak transport substancji był jedyną funkcją węzła, można by 
było uznać węzeł za strukturę nieadekwatną do spełnianej funkcji, ponie- 
waż o wiele sprawniejszy byłby transport w plesze nie podzielonej na ko- 
mórki. Poza rolą transportera węzeł spełnia co najmniej jeszcze jedną rolę, 
nie docenianą w literaturze. Przy uszkodzeniu któregoś z miedzywezli częś- 
ci rośliny oddzielone od uszkodzonej komórki węzłami pozostają żywe. 
Węzeł jest więc strukturą, która w przypadku uszkodzenia jednego z mię- 
dzywezli zabezpiecza sąsiednie komórki przed śmiercią. Ta funkcja węzła, 
jak się wydaje, musiała mieć istotne znaczenie w ewolucji ramienic. 

Uszkodzenie jednego z międzywęźli stwarza niebezpieczeństwo śmierci 
sąsiednich komórek i prawdopodobnie zapoczątkowuje skomplikowane pro- 
cesy reparacyjne w węzłach przylegających do uszkodzonego międzywęźla. 
W pracy rozważany jest tylko problem wpływu objętościowego cytoplazmy. 
Wiadomo, że wewnątrz żywych komórek ramienie panuje pokaźne ciśnie- 
nie hydrostatyczne (turgor) rzędu 6 atmosfer. Z drugiej strony, jak wyka- 
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zaty badania mikroskopowe Spanswicka i Costertona [22] oraz 
Spiewli [23], w ściankach komórek węzła znajdują sie plazmodesmy 
o średnicy 50—60 nm, dzięki którym cytoplazma całej rośliny jest ciągła. 
Jeżeli w pewnym miedzywezlu na skutek uszkodzenia ciśnienie hydrosta- 
tyczne spadnie do poziomu ciśnienia zewnętrznego, to na sąsiednich wę- 
złach wytworzy się różnica ciśnień rzędu 6 atmosfer, co może spowodować 
wypływ objętościowy cytoplazmy poprzez plazmodesmy. Praca jest po- 
święcona zbadaniu możliwości takiego wypływu i ewentualnych mecha- 
nizmów zmniejszających ten wypływ. 


MATEMATYCZNE SFORMUŁOWANIE PROBLEMU 


Kilka komórek międzywęzłowych przedzielonych węzłami przedsta- 
wiono schematycznie na ryc. 1. Załóżmy, że komórka b została uszko- 
dzona (odcięta). Ciśnienie w komórce a bezpośrednio po obcięciu komórki 
b nie ulega zmianie, natomiast w komórce b maleje do zera (dla prze- 
pływu cieczy istotna w tym przypadku jest różnica ciśnień, zatem można 
przyjąć stałe ciśnienie zewnętrzne za równe zeru). Pod wpływem pow- 
stałej różnicy ciśnień p zawartość komórki a wypływa poprzez plazmo- 
desmy węzła na zewnątrz. 


Ryc. 1. Odcinek rośliny Nitellopsis 
obtusa z węzłami: a, b — między- 
węźla, W — węzeł, C — miejsce 
przecięcia międzywęźla b, Vy — 
objętość wody wpływającej do ko- 
mórki a przez powierzchnię bocz- 
ną, Ve — objętość cytoplazmy wy- 
pływającej z komórki a przez plaz- 
modesmy 


Wypływowi temu towarzyszy obniżanie się ciśnienia wewnątrz ko- 
mórki. Na skutek obniżania się ciśnienia następuje kurczenie się ścianki 
komórkowej i osmotyczne ssanie wody ze środowiska do wnętrza komórki. 
Dopływ wody ze środowiska zmniejsza stężenie osmotycznie czynnych sub- 
stancji wewnątrz komórki i powoduje obniżenie ciśnienia osmotycznego 
wewnątrz komórki. Procesy te można opisać układem równań (1—6). 


dV;=kS(II— p)dt (1) 
dV,= Apdt (2) 
dV =dV,—dV, (3) 


Il =cRT (4) 


14 


Wypływ objętościowy cytoplazmy przez plazmodesmy węzłów glonów... 301 


dp= mA (5) 
dc _ _ C dV, 
GR LSK. dł (6) 


gdzie: V, — objętość wody wpływającej ze środowiska do międzywęźla a 
(patrz ryc. 1), Va — objętość cytoplazmy wypływającej z komórki a przez 
węzeł, V — objętość komórki międzywęzłowej a, p — ciśnienie wewnątrz 
komórki a (dokładniej różnica ciśnień między wnętrzem komórki a i ośrod- 
kiem, // — ciśnienie osmotyczne wewnątrz komórki a (dokładniej — różn:- 
ca ciśnień osmotycznych), c — stężenie molowe cząsteczek decydujących o 
wielkości ciśnienia osmotycznego w komórce a, t — czas liczony od chwili 
przecięcia komćrki b, k — współczynnik przepuszczalności wody przez bło- 
nę komórkową, S — powierzchnia komórki międzywęzłowej a, « — moduł 
sprężystości ścianki komórkowej, R — stała gazowa, T — temperatura ab- 
solutna. 

Przeanalizujmy sens fizyczny poszczególnych równań. 

Równanie (1) jest równoważne stwierdzeniu, że strumień wody wpły- 
wającej do komórki jest wprost proporcjonalny do różnicy między ciśnie- 
niem osmotycznym roztworu wypełniającego komórkę a ciśnieniem hydro-- 
statycznym wewnątrz komórki. Współczynnik przepuszczalności k, jak 
wynika z pracy Tazawy i Kiyosawy [25], nie zależy od ciśnienia, 
nie zależy również od tego, czy przepływ wody następuje pod wpływem 
różnicy ciśnień osmotycznych, czy hydrostatycznych. Na temat zależności 
współczynnika k od kierunku przepływu wody zdania są podzielone [2, 11, 
24, 25]. Jednak w rozważanym tutaj przypadku nie ma to znaczenia, gdyż 
przepływ wody następuje tylko w jednym kierunku — do wnętrza ko- 
mórki. Wartość współczynnika k wynosi, zgodnie z danymi zawartymi w 
piśmiennictwie: około 1,089 - 10-11 cm?s/g, [2, 11], 1,347 - 10-11 cm?s/g [25], 
2,11 - 10-11 cm?s/g [24]. 

Równanie (2) wyraża proporcjonalność prędkości przepływu cytopla- 
zmy poprzez plazmodesmy do różnicy ciśnień na węźle. Jest ona konse- 
kwencją założenia newtonowskiego charakteru wypływu. W pracach K a- 
miyi i Kurody [10] oraz Donaldsona [8] stwierdzono, że cyto- 
plazma w komórce międzywęzłowej jest cieczą nienewtonowską i jej prze- 
pływ można opisać modelem potęgowym Ostwalda — de Waele'a (7) lub 
modelem pseudoplastyku Binghama (8) 


pa raj (7) 


Pann gy: T>Ty (8) 
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gdzie r — strumień pędu (siła tarcia wewnętrznego działająca na po- 
wierzchnię jednostkową), u — prędkość, y — współrzędna prostopadła do 
prędkości. 

Do obliczeń bardziej przydatny jest model potęgowy (7), jednak war- 
tość parametru a w tym modelu nie jest jednoznacznie określona, wg 
Donaldsona jest ona 12 razy większa niż według Kamiyi (parametr a wcho- 
dzi do wzoru na prędkość wypływu w trzeciej potędze, tzn. prędkość wy- 
pływu przy zastosowaniu tego modelu byłaby określona z dokładnością 
do trzech rzędów wielkości). Poza tym model potęgowy jest — podob- 
nie jak i inne modele nienewtonowskie — modelem empirycznym i sto- 
suje się tylko w warunkach, w których zostały dopasowane parametry. 
Warunki przepływu cytoplazmy przez plazmodesmy są jednak zupełnie 
różne od warunków ruchu rotacyjnego wewnątrz miedzywezla, do których 
dopasowano model. W tej sytuacji zdecydowałem się na przyjęcie założe- 
nia, że cytoplazma w plazmodesmach zachowuje się jak ciecz newto- 
nowska. 

Przy powyższym założeniu stała A wyraża się wzorem (9) 


nrt 
8 nL (9) 


A=N 


gdzie: N — liczba plazmodesm, r — średni promień plazmodesmy, 7 — 
współczynnik lepkości cytoplazmy, L — długość plazmodesmy. Wzór (9) 
jest bezpośrednią konsekwencją prawa  Hagena-Poiseuille'a. War- 
tość stałej A oszacowano na podstawie danych mikroskopowych [22, 23], 
współczynnik lepkości cytoplazmy 70,05 g/cm*s wzięto z pracy P e- 
karka [17]. Przyjęty współczynnik lepkości jest charakterystyczny dla 
ciekłej fazy cytoplazmy (bez inkluzji). Zatem jest to praktycznie dolna 
granica dla wartości współczynnika lepkości. Należy więc liczyć się z 
możliwością, że obliczona w oparciu o taką wartość współczynnika lepkości 
stała A jest zawyżona. Również traktowanie plazmodesm jako walcowa- 
tych otworów prowadzi do zawyżenia stałej A. Wypływ obliczony na 
podstawie równań (1—6) będzie więc najprawdopodobniej większy od 
rzeczywistego. Wartość stałej A oszacowano na 1,12: 10711 cm?s/g. 

Równanie (2) opisuje przepływ stacjonarny. Pominięcie procesów 
przejściowych nie wnosi praktycznie żadnego błędu do wyniku analiz, po- 
nieważ charakterystyczny czas tych procesów ze względu na małe roz- 
miary plazmodesm jest rzędu 107% sek. [3]. Równania (3 i 4) nie wyma- 
gają objaśnień. 

Równanie (5) jest „równaniem stanu” ścianki komórkowej i stwier- 
dza, że zmiana ciśnienia w komórce jest proporcjonalna do względnej 
zmiany objętości. Liniowość tej zależności jest w badanym przypadku 
dobrym przybliżeniem, gdyż chodzi tu o małe zmiany objętości (maksi- 
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mum 1%). Wartość modułu sprężystości waha się, według różnych auto- 
rów [2, 24], w granicach 10%—10% dyn/cm?. Ostatnie równanie (6) opisuje 
zachodzącą podczas wypływu objętościowego zmianę stężenia roztworu 
wewnątrzkomórkowego. Postać matematyczna tego równania nie jest tak 
oczywista jak poprzednich (1—5), dlatego też niżej podane jest jego 
wyprowadzenie. Zmianę stężenia w ciągu skończonego dostatecznie ma- 
łego czasu At można wyrazić następująco: 


dow RÓB. 
V+4V,—4V> 
czyli 
AV 

dex — ZE 

rd V+4V,—4AV, 
stąd 

Ac c AV, 


At  V+AVY SMY, at 
Przechodząc w ostatnim wzorze do granicy przy At > 0 otrzymujemy 


de c dV; 


dt V dt (6) 


ROZWIĄZANIE UKŁADU RÓWNAŃ (1—6) 


1. Rozwiązanie numeryczne 


Układ równań (1—6) po przekształceniach sprowadza się do równania 
drugiego rzędu dla jednej zmiennej — ciśnienia (patrz przypis): 


dp 


ze | + fp’ = (10) 


(p+) SR + (yp-+90) SB + | 
e 
y= y, 8A+2k5) 
= —+_(A+ks) 
+ Vo 
gia 
p= Se A(A+kS) 


ĝ= €&—Po 


Vo, Po— objętość komórki a i ciśnienie wewnątrz komórki a przed prze- 
cięciem komórki b (ryc. 1). Dla Vo i po przyjęto wartości: 
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Vo=9,815 * 1073%cm3 
Po=6 - 105 dyn/cm? 


Mamy wiec 
Ple- o= Po (11) 
Poza tym należy zbadać wartość SP ala chwili t=0. Z równania (5) wy- 
nika, że 
pE TATOV 
dt Vo dt (12) 
Wstawiając (3) do (12) otrzymamy: 
dt zł dti dt | ag 
i wykorzystując równania (1) i (2) 
- SRE Acha 
Fi Cv, [kS (II — p) — Ap] (14) . 
Dla chwili t=0 zachodzi /J=p=pp, zatem 
dp eAp, 
-= - = = 1 
Capea Vo (LĄ 


Równanie (10) wraz z warunkami początkowymi (11) i (15) jedno- 
znacznie określa funkcję p(t). Niestety jest ono równaniem nieliniowym 
i jego rozwiązania nie można otrzymać w postaci analitycznej. Dlatego 
też po uprzednim sprowadzeniu tego równania do układu dwóch równań 
pierwszego rzędu, zostało ono rozwiązane numerycznie. 

Przy określonych wartościach pozostałych parametrów (ryc. 2) ciśnie- 
nie wewnątrz komórki po obcięciu sąsiedniego międzywęźla spada w cią- 
gu 5—10 sek do wartości 2,5 * 10% dyn/cm?, co stanowi około 40% wartości 
początkowej i maleje dalej bardzo wolno. Od wartości e zależy tylko szyb- 
ko malejąca część rozwiązania. Wartość modułu sprężystości określa więc 
czas relaksacji ścianki komórkowej. Czas ten jest stosunkowo krótki i pro- 
cesy zachodzące po upływie 10 sek. nie zależą od własności sprężystych 
ścianki komórkowej. 

Zależność rozwiązania p(t) od wartości współczynnika przepuszczal- 
ności k ma nieco inny charakter (ryc. 3) i jest ono niemal jednakowe 
dla obu wartości k w ciągu pierwszych sekund. O zmianach ciśnienia w 
komórce w ciągu pierwszych kilku sekund po obcięciu sąsiedniej ko- 
morki, decydują tu własności sprężyste ścianki (ryc. 2), natomiast wpływ 
zjawisk osmotycznych w tej fazie jest znikomy (ryc. 3). Z drugiej 
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p [10° dyn/cm] 


Po" 6:10Śdyn/em? 
5 fe he 
ZM crą?5 
k= 1089910) = 
P £ V,=9,815« 1073%cm3 
2,5110 


S=0,785cm? 


10 20 30 40 50 60 70 80 . 90 


Ryc. 2. Zależność funkcji p (t) od wartości parametru e 


P [10 dyn/cm] 

6 
E= 5110? dyn/cm? 

5 B= 610° dyn/cm? 

-11 emes 

A=1,12:10 SSG 

4 V;=9815:10%cm* 
S=0,785 cm? 


k=1bx10 " 


k=1,08910" 


10 20 30 40 50 60 70 80 90 


Ryc. 3. Zależność funkcji p (t) od wartości parametru k 


100 t [s] 


100 t [5] 
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strony po upływie dostatecznie długiego czasu wartość ciśnienia wewnątrz 
komórki przestaje zależeć od modułu sprężystości ścianki (ryc. 2), ale wy- 
raźnie zależy od prędkości osmotycznego ssania wody (ryc. 3). Fakt ten 


20 — Annales 
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daje możliwość uzyskania przyblizonych rozwiązań analitycznych układu 
równań (1—6). 


2. Przybliżone rozwiązanie analityczne 
dla fazy początkowej 


Przybliżony opis początkowej fazy procesu wypływu cytoplazmy moż- 
na uzyskać pomijając całkowicie zjawisko osmozy, tj. wprowadzając do 
układu równań (1—6) dodatkowy warunek 


k=0 (16) 
Wówczas układ równań (1—6) zredukuje się do równania (17): 
dp „1YSA 
a av? (17) 


z warunkiem początkowym (11), którego rozwiązaniem jest funkcja: 


tA 


p=poe ve (18) 


Z (18) w oparciu o równanie (2) można wyznaczyć objętość cytoplazmy 
wypływającej przez plazmodesmy: 


vale Let.) | (19) 


3. Przybliżone rozwiązanie analityczne 
po upływie dostatecznie długiego czasu 


Rozwiązanie równania (10) zależy od wartości modułu sprężystości 
ścianki komórkowej tylko przez parę pierwszych sekund (ryc. 2). W tym 
czasie następuje wypływ rzędu 107% cm’. Jest to objętość bardzo mała w 
porównaniu z objętością komórki (9,815 * 107% cm?), ale wystarczająco du- 
ża dla zrelaksowania ścianki komórkowej. O dalszym przebiegu procesów 
decyduje stosunek strumieni: wpływającego przez powierzchnię boczną 
i wypływającego przez plazmodesmy. A zatem, jak się wydaje, można 
przyjąć założenie, że objętość komórki pozostaje stała i że proces sprę- 
żystej relaksacji ścianki komórkowej jest natychmiastowy. Przy tych za- 
łożeniach dV=0, a układ równań (1—6) upraszcza się do następującej 
postaci: 

dV;=kS(II —p)dt (20) 


dV;= Apdt (21) 
dV,=dV, (22) 
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II=cRT : (23) 
JEŻY GydVi 
dt, V» dt 183 


gdzie V, jest stałą objętością komórki. 
Wstawienie (20) i (21) do (22) daje równanie algebraiczne: 


kS(II— p)= Ap (25) 
a wstawienie (20) i (23) do (24) daje równanie różniczkowe: 
dll H 
E y,*SUT- p) (26) 
Z (25) wynika, że 
kS 
P=KS+A (27) 


Po wstawieniu (27) do (26) orzymujemy równanie różniczkowe z jedną 
zmienną: 


SA | AD TRSA 

dt V, kS+A 
Jest to równanie nieliniowe, ale łatwo jest uzyskać jego rozwiązanie ana- 
lityczne drogą rozdzielenia zmiennych. Formułując warunek początkowy, 
nie można w tym przypadku przyjmować, że dla t=0 ciśnienie wewnątrz 
komórki jest równe ciśnieniu przed obcięciem sąsiedniej komórki (py), 
gdyż przeczyłoby to założeniu o natychmiastowej relaksacji ścianki. Wiel- 
kością, która w początkowej chwili wypływu będzie miała taką samą war- 
tość jak przed obcięciem sąsiedniego międzywęźla, jest ciśnienie osmo- 
tyczne wewnątrz komórki. Zatem warunek początkowy będzie miał po- 
stać: 


(28) 


IT, 4F7Po (29) 
Równanie (28) można przekształcić do następującej postaci: 
ot adt (30) 
gdzie 
uae PARA 
Vo (kS+ A) 
Scalkowanie (30) daje zależność: 
— = —at+B (31) 


gdzie B jest stałą całkowania. 
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Wstawienie warunku początkowego (29) do (31) daje 


B=——. 
Po 
Zatem 
groet Basle, 
Pott+1 (32) 
Zależność p(t) można otrzymać wstawiając (32) do (27) 


i Bogny Pr 
PO kS+A patti 
Objętość cieczy wypływającej przez plazmodesmy (V>) otrzymać moż- 


‘na wstawiając wartość p (33) do wzoru (34), który jest konsekwencją 
równania (21): 


(33) 


=A fp()dr (34) 
0 d 


Z przekształceń tych otrzymujemy: 
Va =Vo In (1 + poat) (35) 
Wzory (32), (33) i (35) stanowią uproszczone rozwiązanie modelu. 
Z porównania wynika, że funkcja (18) jest dobrym przybliżeniem 
rozwiązania dla t < 0,5s, natomiast funkcja (33) jest dobrym przyblize- 
niem rozwiązania dla t > 6s (ryc. 4). 


P [Mo*dyn/em?] - 

6 
£=5"40°dyn/em? 
R610 dyn/cm? 

5 ! Ast saxon ERES 
k *1,432:10 tents 

4 V*9,815210 Sem? 

j S =0785 cm? 


10 20 30 40 50 60 70 80 90 100 t [s] 


Ryc. 4. Rozwiązanie układu równań (1—6): 1 — numeryczne, 2 — przybliżone dla 
fazy początkowej, 3 — przybliżone dla dostatecznie dużych czasów l 
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DYSKUSJA 


1. Kryterium czasu 


Jeżeli lepkość cytoplazmy miałaby zabezpieczyć komórkę a (ryc. 1) 
przed wypływem objętościowym po odcięciu komórki b, to bardzo ważne 
jest określenie czasu, w jakim ma trwać to zabezpieczenie. Zarówno z na- 
tury ruchu cieczy lepkiej, jak i z obu rozwiązań modelu matematycznego 
wynika, że przy otwartych plazmodesmach wypływ objętościowy może 
być dowolnie duży, jeżeli będzie trwał dostatecznie długo. 

W nienaruszonej roślinie cytoplazma nie kontaktuje się ze środowis- 
kiem bezpośrednio, lecz jest od niego oddzielona plazmalemmą i ścianką 
komórkową. Dla zabezpieczenia przed wypływem objętościowym bardziej 
istotną strukturą jest ścianka komórkowa, której sprężystość równoważy 
panujące wewnątrz komórki ciśnienie. Zachodzące po uszkodzeniu rośliny 
(przecięciu międzywęźla) procesy reparacyjne powinny prowadzić do przy- 
wrócenia naturalnej sytuacji, w której protoplast jest oddzielony od śro- 
dowiska ścianką komórkcwv49, tj. do zabudowania światła plazmodesm sub- 
stancją ścianki. Substancją .zcydującą o dużej wytrzymałości mechanicz- 
nej ścianki jest celuloza. Włókna celulozy wytrzymują napięcia do 
110 kG/mm? [18]. Na podstawie tej wartości i rozmiarów plazmodesm 
można obliczyć, że dla zrównoważenia panującego w komórce ciśnienia 
wystarczy warstwa celulozy o grubości kilku A, a zatem wystarczająca 
jest jedna warstwa celulozy krystalicznej. Opierając się na rozmiarach 
elementarnej komórki celulozy krystalicznej (10,3X7,9X8, 35 A) [18] 
można obliczyć, że dla przesłonięcia światła plazmodesmy jedną warstwą 
celulozy krystalicznej (najgęściej upakowanej) powinny być spolimery- 
zowane łańcuchy celulozy o łącznej długości 105 A (2-104 cząsteczek glu- 
kozy). 

Jednak literatura dotycząca ścian komórkowych nie zawiera bliższych 
informacji o syntezie celulozy [4, 15, 19]. Informacją, którą można wyko- 
rzystać do oszacowania czasu zasklepienia plazmodesm celulozową ścian- 
ką, są wyniki badań Millmana i Colvina [14], dotyczące syntezy 
celulozy na powierzchni agaru przez Acetobacter xylinum. W pracy tej 
stwierdzono między innymi, że mikrofibryle celulozy wydłużają się z 
prędkością 420 A/s. Jeżeli przyjąć, że na każdą plazmodesmę przypada 
jedna aktywna cząsteczka enzymu, to zasklepianie plazmodesmy powinno 
trwać 4—5 min. Czas ten może być w rzeczywistości nieco krótszy, jeżeli 
stężenie enzymu jest większe. Może jednak być o wiele dłuższy, gdyż nie 
jest wykluczone, że w tym przypadku o szybkości reakcji decyduje dy- 
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fuzja, a nie aktywność enzymu. Tym niemniej powyższe oszacowanie daje 
orientację co do rzędu wielkości czasu. 

Innym możliwym procesem reparacyjnym jest odkształcenie ścianki 
na skutek zmniejszenia się — po odcięciu międzywęźla — ciśnienia hy- 
drostatycznego w plazmodesmach kontaktujących się ze środowiskiem. 
Gdyby odkształcenie ścianki było dostatecznie duże, mogłoby ono prowa- 
dzić do zamykania plazmodesm. Jednak natychmiastowe (sprężyste) de- 
formacje ścianki u Nitelli nie przekraczają 5% [5, 12]. Jest to odkształ- 
cenie zbyt małe. Natomiast dalsze odkształcenie, powodowane zmianą 
ciśnienia, nie jest procesem czysto fizycznym, lecz płynięciem wiskoela- 
stycznym, któremu towarzyszą reakcje chemiczne (tzw. chemical creep) 
[5]. Ścianka komórkowa wykazuje zmianę własności mechanicznych pod 
wpływem bodźców chemicznych po upływie czasu od kilku do kilkunastu 
minut [5, 18] od chwili zadziałania bodźca. 

Obserwacje zmian oporu elektrycznego węzła po uszkodzeniu rośliny 
[21] pozwalają przypuszczać, że plazmodesmy zostają przesłonięte błoną 
po 10—25 min. od obcięcia międzywęźla. Wydaje się więc, że w ciągu 
tego czasu procesy reparacyjne doprowadzają co) całkowitego zatrzyma- 
nia wypływu objętościowego cytoplazmy, co ¿aje możliwość wytworze- 
nia na powierzchni cytoplazmy błony selektywnej. 

W dalszych rozważaniach nad wynikami modelu matematycznego za- 
kłada się, że komórka jest zabezpieczona przed wypływem objętościowym, 
jeżeli w ciągu 10 min. od obcięcia sąsiedniej komórki zmiany powodowa- 
ne wypływem objętościowym są dostatecznie małe. 


2. Plazmodesmy jako otwarte kanały 


Załóżmy, że plazmodesmy są otwartymi kanałami o wymiarach otrzy- 
manych w badaniach mikroskopowych [22, 23]. Przyjmijmy, że komórka a 
na ryc. 1 ma długość 5 cm i średnicę 0,5 mm. Komórka o takich wymia- - 
rach ma powierzchnię s=0,785 cm? i objętość Vą*9,815 * 107% cm*. Inne 
parametry zostały omówione wyżej. Ich wartości są równe: współczynnik 
przepuszczalności błony komórkowej dla wody k=1,432 : 10711 cm? s/g, 
ciśnienie wewnątrz nienaruszonej komórki pọ=6 -10° dyn/cm?, stała cha- 
rakteryzująca prędkość wypływu przez plazmodesmy A= 1,12:10-11 cm‘ s/g. 
Przy tych wartościach parametrów współczynnik a, występujący we wzo- 
rach (32), (33) i (35) jest równy 5,717 : 10-70 cm s/g, pya= 3,43 * 1073%s-1 i dla 
czasu t=10 min=600 s, pyat=2,058. Ze wzorów (32, 33, 35) wynika, że 
w ciągu 10 min. ciśnienie osmotyczne wnętrza komórki spadałoby do 1/3 
początkowej wartości, ciśnienie hydrostatyczne do 16,5% wartości począt- 
kowej, a przez plazmodesmy wypłynęłaby objętość równa 112% objętości 
miedzywezla. 
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Jakkolwiek nie ma ilościowego kryterium śmierci komórki i trudno 
jest przesądzać, czy powyższe zmiany spowodowały śmierć komórki, 
pewne jest jednak, że takie zmiany nie zachodzą w rzeczywistości, bo- 
wiem komórki międzywęzłowe ramienic po odizolowaniu od pozostałej 
części rośliny zachowują turgor i wykazują normalny ruch cytoplazmy. 
Potraktowanie plazmodesm jako otwartych kanałów prowadzi zatem do 
sprzeczności z doświadczeniem. Podobna sprzeczność występuje — przy 
założeniu otwartych kanałów plazmodesm — między wartością oporu 
elektrycznego węzła otrzymaną doświadczalnie a wartością teoretyczną 
(patrz niżej). 


3. Hipoteza o zablokowaniu częściowym 
plazmodesm 


Załóżmy, że przekrój plazmodesmy wygląda w sposób przedstawiony 
na ryc. 5, tzn. centralną część plazmodesmy zajmuje nieruchoma struk- 
tura, cytoplazma zaś wypełnia pozostającą, wolną część plazmodesmy. 
Niech R, oznacza opór elektryczny plazmodesmy całkowicie wypełnionej 
cytoplazma, a Rə opór plazmodesmy przedstawionej na ryc. 5. 


Ryc. 5. Przekrój częściowo zablokowanej plazmodesmy: r — 
promień plazmodesmy, x r — promień umieszczonej cen- 
tralnie struktury blokujacej; 0<~ <1 


Wtedy 


ŻA eh nig (36) 


Jeżeli oznaczymy przez Q; objętość cieczy wypływającej w ciągu 1 sek. 
przez plazmodesmę otwartą, a przez Q, — plazmodesmę częściowo zablo- 
kowaną, to: 


Q: (1 — x*) 
— =(1—x) [3] (37) 
Qı In A 
x 
Poniewaz 
Q: _ A> 
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gdzie: A, — stała A (2) charakteryzująca wypływ przez węzeł, którego 
plazmodesmy są otwarte, A> — analogiczna stała dla węzła, którego plaz- 
modesmy są częściowo zablokowane, to otrzymamy: 


A> (1 — x?) 
4=(1—x)— —- (38) 
tt In 


Zależność stosunków Ry/Rą i A/A, od x jest przedstawiona na ryc. 6. W 
związku z tym, że dla x > 0,9 stosunek A;/A, jest już rzędu błędu popeł- 
nianego przy obliczeniach z pomocą tablic czterocyfrowych, na ryc. 7 
przedstawiono interpolację liniową obu funkcji między x=0,9 i x=1. 
Dla x=0,9, R;/Ra=0,19, A,/A,;=0,0013. Dla x=1 obie funkcje mają war- 
tość 0. 
R4/R2 
A>/A4 : 

10 


R, /R2 


01 Ryc. 6. Zależność oporu plazmodesmy 
i prędkości przepływu od promienia 
01 02 03 0% OS 06 07 08 09 10x struktury blokującej 


Jak wynika z badań innych autorów [22, 23] oraz z pomiarów prze- 
prowadzonych w naszej pracowni [20], rzeczywisty opór plazmodesm jest 
od kilkudziesięciu do kilkuset razy większy od obliczonego na podstawie 
danych mikroskopowych. Nie będzie więc przesadą założenie, że Ry/R+=0,1. 
Wtedy, jak wynika z ryc. 7, As/A;=0,00069, czyli A»=4; * 6,9 * 10-4. Za- 
tem założenie, że częściowe zablokowanie plazmodesm zwiększa ich opór 10 
razy, oznacza, że wypływ objętościowy powinien być o trzy rzędy wielkości 


mniejszy, a stała A powinna być równa 7,728: 10% cm‘s/g. Przy tych 


samych wartościach pozostałych parametrów co poprzednio współczyn- 
nik a=7,866 : 10-13 pya=4,72 *1075s71. Dla czasu t=10 min = 600s, 
poat =2,83 - 10-3. Wstawiając tę wartość do wzorów (32, 33, 35) otrzy- 
mamy: 


= um m 
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Ryc. 7. Interpolacje liniowe R,/R, i A,/A, 
w przedziale 0,9 Sx S1 


IT 1 


p kS 


Po  1+poat 
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R,/tS 


020 


018 


016 


0,14 


012 


010 


qos 


006 


0,04 


0,02 


= 0,997 


1 


Pp kS+A 1+p,at 


Və 


0 


= 0,996 


= In (1 +poat) = 0,003 


Powyższe wartości oznaczają, że w ciągu 10 min. zachodzi zmniejszenie 
ciśnienia osmotycznego o 0,3% w stosunku do wartości początkowej 
i zmniejszenie ciśnienia hydrostatycznego o 0,4% oraz wypływ poprzez 
plazmodesmy objętości równej 0,3% objętości komórki międzywęzłowej. 
Powyższe zmiany zachodzące w międzywęźlu po obcięciu sąsiedniego 
miedzywezla są bardzo małe. O tym, że nie stanowią one zagrożenia dla 
życia komórki świadczy fakt, że tego samego rzędu zmiana ciśnienia 
i wypływ objętościowy towarzyszą każdemu potencjałowi czynnościowe- 
mu, generowanemu przez komórkę [1, 2]. 
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WNIOSKI KONCOWE 


Przedstawione wyniki modelu matematycznego nie odpowiadają do- 
kładnie rzeczywistym zmianom zachodzącym w komórce a po odcięciu 
komórki b (ryc. 1). Przyczyną tego jest między innymi charakteryzowa- 
nie wielkości wypływu przez węzeł stałym parametrem A. Jakakolwiek 
byłaby natura procesów reparacyjnych, powinny one prowadzić do 
zmniejszenia wypływu objętościowego. Zatem w modelu matematycz- 
nym parametr A powinien być nie wielkością stałą, a malejącą funkcją 
czasu. W rzeczywistości więc Vs powinno być mniejsze od otrzymanych 
wartości. Mniejsze powinny być również spadki ciśnienia osmotycznego 
(x) i turgoru (p). Jeżeli więc z rozwiązania równań (1—6) wynika, ze 
zmiany zachodzące w komórce są na tyle małe, że nie zagrażają jej życiu, 
to można uważać, że zabezpieczenie jest dostateczne, ponieważ rzeczy- 
wiste zmiany są jeszcze mniejsze. 

Z drugiej strony, jak wynika z własności ścianki komórkowej [5, 13, 
16], nie istnieją mechanizmy reparacyjne o natychmiastowym działaniu. 
Zatem lepki charakter wypływu w połączeniu ze strukturą węzła wy- 
daje się niezbędnym elementem w procesach prowadzących do zachowa- 
nia. komórki. Należy jednak zaznaczyć, że mechanizm wypływu lepkiego 
jest zabezpieczeniem dostatecznym tylko przy przyjęciu hipotezy o częś- 
ciowo zablokowanych plazmodesmach. 

Hipoteza o częściowym zablokowaniu plazmodesm opiera się w głów- 
nej mierze na wynikach pomiarów oporu elektrycznego węzła. Ani po- 
miary elektryczne, ani przeprowadzone rozważania nie przesądzają jed- 
nak o tym, w jaki sposób plazmodesmy są zablokowane. Możliwe, że 


plazmodesmy są blokowane w inny sposób aniżeli przedstawiono to na ~ 


ryc. 5 (np. przez organelle komórkowe lub inne stosunkowo sztywne in- 
kluzje cytoplazmy). Bez względu jednak na sposób blokowania przyjęcie 
hipotezy o częściowym zablokowaniu wyjaśnia rozbieżności między wy- 
nikami pomiarów elektrycznych i mikroskopowych obserwacji i jedno- 
cześnie pozwala zrozumieć działanie mechanizmu zabezpieczającego ko- 
mórkę przed wypływem objętościowym. 


PRZYPIS 


Sprowadzenie układu równań (1—6) do równania z jedną zmienną. 
Równania (1) i (2) wstawiamy do (3) i otrzymujemy: 


dv 


gt = KSI- (ks +A)p_ (39) 
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Z równania (6) i (1) wynika (40): 


dy e 
dt ~~ yKSUI-p) (40) 
Do (39) i (40) wstawiamy zależność II =cRT i otrzymujemy (41) i (42): 
dv 
aor kScRT — (kS+4A)p (41) 
de e 
dt 7 y KS(ERT—p) (42) 
Z równania (5) wynikają zależności: 
dV _ Vo dp 
dt 7 a dt (43) 
Vo 
V=- (P—Po) + Vo (44) 


Po wstawieniu (43) i (44) do (41) i (42) otrzymujemy układ dwóch równań róż- 
niczkowych ze zmiennymi pic 


dp e e ; 
RA kScRT — y (kS+A)p (45) 
dc ec 
“Gt vipe pj tS (RT — p) (46) 
Z (45) po rozwikłaniu otrzymujemy zależnyść c od p i s ` 


kS+A Vo dp 
C= ESRT ?* skSRT dt 
Skąd po zróżniczkowaniu otrzymujemy: 
dc kS+A dp Vo dp 
dt kSRT dt ' ekSRT dt? (48) 
Po wstawieniu (47) i (48) do równania (46) i wykonaniu przekształceń alge- 
braicznych otrzymujemy ANO | 


(47) 


dp 
+0 gg + ptg) Gr + (2 2y +8p’=0 (10) 
gdzie: ô = £ — Po 


p= y (GA +2k8) 
g=q, (A+kS) 


at 
e-(y, ] A(A+KkS) 
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PESIOME 


B pa6ore cdbopMyJrupoBaHa u pellieHa cucremMa „rucbcbepeHquaJIbHBIX 
ypaBHeHuii, ONUCIBAIOLIMX OObEMHOe BbITeKAHMe IUTONJIAZMbI IPA JIO- 
KAJIBHOM IIOBpezK/JeHMM paCTeHuA Characeae. Pe3yJIbTaTbI pacH6TOB yKa- 
3BIBAIOT Ha TO, YTO IJIa3MOJECMBI B y3JAx BojqopocJeń Characeae sa- 
CTMUHO OJIOKMPOBAHBI. 


SUMMARY 


A system of differential equations has been formulated and solved. 
The system describes the volume outflow of cytoplasm during a local 
injury of the Characea plant. The results of the calculations suggest that 
plasmodesmata in the nodes of Characeae algae are partly blocked up. 


